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PREFACE 


The  design  of  hypersonic  vehicles  has  stimulated  a  large  number 
of  fluid-mechanical  studies  concerned  with  the  hypersonic  boundary 
layer.  This  Memorandum  enlarges  the  scope  of  these  studies  to  include 
various  models  for  the  fluid  viscosity,  the  effect  of  high  mass-injec¬ 
tion  rates  at  the  stagnation  point,  and  solutions  at  the  hypersonic 
limit  where  the  free-stream  velocity  is  very  high.  These  results 
should  be  useful  as  a  guide  in  the  design  of  vehicles  operating  at 
hypersonic  speed. 

This  study  is  one  of  a  series  dealing  with  the  hypersonic  laminar 
boundary  layer,  and  part  of  continuing  RAND  research  in  theoretical 
fluid  mechanics. 
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SUMMARY 


This  Memorandum  presents  certain  numerical  solutions  to  the  lami¬ 
nar  boundary- layer  similarity  equations  which  illustrate  the  effects 
of  Prandtl  number,  viscosity- temperature  variation,  mass  transfer, 
wall  temperature,  pressure  gradient,  and  hypersonic  parameter  a  on  the 
structure  of  the  boundary  layer  and  the  derivatives  of  velocity  (shear 
force)  and  total  enthalpy  at  the  wall. 

The  results  show  that  the  common  boundary- layer  simplifications 
(such  as  Prandtl  number  *  1,  viscosity  proportional  to  temperature, 
zero  pressure  gradient)  often  lead  to  numerical  errors  of  20  to  50  per 
cent  in  the  predicted  heat  transfer,  skin  friction,  and  displacement 
thickness.  The  shear  and  enthalpy  gradient  are  shown  to  depend  on  the 
form  of  the  viscosity- temperature  relation,  the  Prandtl  number,  and  the 
hypersonic  parameter  a.  Results  for  high  surface  mass -injection  rates 
at  the  stagnation  point  indicate  that  the  shear  at  the  wall  does  not 
vanish  as  in  the  case  of  a  flat  plate. 
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I.  INTRODUCTION 


In  most  physical  situations,  the  laminar  boundary  layer  which  de¬ 
velops  along  the  surface  of  a  body  in  high-speed  flow  does  not  behave 
in  an  exactly  similar  manner.  That  is,  the  mathematical  restrictions 
necessary  to  reduce  the  partial  differential  equations  governing  bound¬ 
ary-layer  motion  to  ordinary  differential  equations  in  transformed  co¬ 
ordinates  are  not  exactly  satisfied.  Of  the  several  concepts  which 
have  been  proposed  to  deal  with  such  circumstances,  local  similarity 
is  perhaps  the  most  successful.  The  local  similarity  concept  assumes 
that  at  every  streamwise  station  the  boundary  layer  adjusts  to  changes 
in  the  geometric  and  thermodynamic  boundary  conditions  and  is  identical 
in  all  essential  respects  to  the  similar- solution  boundary  layer  whose 
history  includes  the  local  boundary  conditions. 

In  this  context,  similar  solutions  to  the  laminar  boundary- layer 
equations  are  of  great  utility  in  predicting  the  skin  friction,  heat 
transfer,  and  boundary- layer  displacement  thickness  for  a  wide  variety 
of  "smooth*1  body  shapes.  A  large  number  of  solutions  to  the  low-speed 
and  stagnation-point  boundary- layer  problems  have  been  published,  but 
little  systematic  work  has  been  done  on  the  important  effects  of  vis¬ 
cous  dissipation  with  Pr  +  1  and  variations  of  the  (pn)  product  with 
temperature.  The  lack  of  adequate  solutions  is  particularly  evident 
in  the  hypersonic  limit,  a  =  (u^/2He)  (u^/u^)^  -*  1. 

This  Memorandum  presents  selected  numerical  solutions  to  the  lam¬ 
inar  boundary- layer  similarity  equations  which  illustrate  systemati¬ 
cally  the  effects  of  Prandtl  number,  viscosity- temperature  variation, 
mass  transfer,  wall  temperature,  pressure  gradient,  and  hypersonic 
parameter  o  on  the  structure  of  the  boundary  layer  and  the  derivatives 
of  velocity  and  total  enthalpy  at  the  wall.  It  will  be  found  that  the 
comnon  boundary- layer  simplifications  (such  as  Pr  *  1,  PM-  »  constant, 
or  3  *  0)  often  lead  to  numerical  errors  of  20  to  50  per  cent  in  the 
predicted  heat  transfer,  skin  friction,  and  displacement  thickness. 


* 

The  quantity  a  will  be  referred  to  as  the  "hypersonic  parameter ,M 
and  o  -*  1  will  be  called  the  "hypersonic  limit. " 


-2- 


Furthermore,  the  shear  and  enthalpy- grad lent  profiles  in  the  boundary 
layer  depend  strongly  on  the  assumed  form  of  the  viscosity- temperature 
relation  and  the  Prandtl  number,  and  also  on  the  hypersonic  parameter 
a . 
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II,  THE  SIMILARITY  EQUATIONS 

The  partial  differential  equations  governing  the  two-dimensional 

compressible  laminar  boundary  layer  in  a  perfect  gas  may  be  reduced  to 

(1)* 

the  following  pair  of  coupled  ordinary  differential  equations: 


+  ff-  -  0  [f,z  ■  (i  •  g  e  •  J  =  o  (i) 


*K'+3riL(1-fc> 


,2 


-  o.  (2) 


a  -  v 


The  boundary  conditions  are  written 


T|  «  0  f'  a  9  ■  0;  f  ■  f  ■  constant  (3) 

T|  -  «  f  -  0  -  1,  (4) 


where  f'  is  the  normalized  velocity  (u/u^) ,  0  is  the  normalized  total 
enthalpy  (H  -  HJ)/ (Hfi  -  H^,  and  o  is  the  hypersonic  parameter 
(u^/2H  )  (ue/uB>)2.  The  derivation  leading  to  Eqs.  (1)  and  (2)  assumes 
that  both  the  specific  heat  cp  and  the  pressure-gradient  parameter  0 
are  constant;  0  is  defined  by  the  equation 


0 


__1  j>,A,u>dx 

V  VwV 


X 


(5) 


★The  effects  of  sweep  may  also  be  studied  in  a  similar  sianner; 
this  increases  the  number  of  coupled  equations  to  three,  and  analogous 
results  may  be  obtained.  The  effects  of  sweep  may  be  found  in  Refs. 

1  and  2. 
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The  similarity  variable  T]  is  defined  by 


TKx.y) 


U 

e 


(2§)1/2 


y 

J  p  dy> 


o 


where 


x 

§(x)  =  J  P^l^dx. 

o 


(6) 


(7) 


Equations  (1)  -  (4)  have  been  solved  using  standard  numeric al- 
integration  techniques.  An  extensive  tabulation  of  solutions  is  pres¬ 
ently  being  prepared  by  the  authors.  The  typical  results  reported  in 
this  Memorandum  were  selected  to  illustrate  the  most  important  conclu¬ 
sions  derived  from  these  solutions.  Empirical  techniques  for  estimat¬ 
ing  effective  fluid-property  variations  across  the  boundary  layer  are 
also  suggested. 

The  physical  significance  of  the  hypersonic  parameter  can  be  seen 
by  noting  its  structure  and  position  in  the  energy  equation.  In  Eq. 
(2),  it  is  clear  that  for  Pr  »  1 ,  the  viscous-dissipation  term  contain¬ 
ing  a  vanishes.  For  Pr  f  1 ,  the  value  of  o  determines  to  some  degree 

the  role  of  the  viscous-dissipation  term  in  the  energy  equation.  The 

2  2 

hypersonic  parameter  a  is  composed  of  two  terms,  (u^/ 2H  )  and  (u  /ua)  . 

2  e  e 

At  the  stagnation  point,  (u  /uj)  *  0;  on  the  other  hand  high-speed 

e  2 

boundary  layers  are  characterized  by  (u  /u  )  ^  1  if  the  free-streaxn 

2  ® 

velocity  is  high,  that  is,  (u^/21^)  -  1. 
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III.  EFFECTS  OF  THE  TEMPERATURE -VISCOSITY  LAW 


Equations  (1)  and  (2)  are  greatly  simplified  if  the  product  (pp.) 
is  constant  across  the  boundary  layer.  Since  p  *  p(x)  only  in  the 

U) 

boundary- layer  approximation,  this  requires  p.  ~  T  with  U)  *  1 .  It 

(3) 

has  been  shown  by  several  authors  (e.g.,  Van  Driest  ')  that  the  form 
of  the  viscosity- temperature  relation  plays  an  important  role  in  de¬ 
termining  the  properties  of  the  boundary  layer. 

Two  alternative  forms  of  p-(T)  have  been  proposed.  These  are  the 
power-law  representation 


uu  1 

p,~T,±£uj£l 

and  the  Sutherland  formula 


(8) 


M- 


(t) 


3/2 


lH-9' 


0  *  s, 


(9) 


where  s  and  uu  are  parametric  constants,  and  t  *  (T/Tq)  . 

In  view  of  the  fact  that  Sutherland's  law  represents  an  approxi¬ 
mation  to  the  intermolecular  potential,  the  power- law  variation  of 
viscosity  appears  to  be  equally  suitable  as  a  one-parameter  represen¬ 
tation  of  the  viscosity- temperature  law,  and  it  offers  certain  con¬ 
ceptual  advantages  when  dealing  with  air  and  other  gases  at  high  tem¬ 
perature.  The  power- law  representation,  however,  provides  a  less 
accurate  representation  of  the  function  n(T)  over  an  extended  range 
of  temperature. 

Two  related  problems  will  be  considered  with  respect  to  the  func¬ 
tion  ti(T)  .  First,  we  wish  to  examine  the  effects  of  the  temperature- 
viscosity  parameters  uu  and  s  on  the  velocity  and  enthalpy  profiles 
within  the  boundary  layer.  And  second,  we  wish  to  demonstrate  that 
the  two  parametric  forms  of  the  viscosity- temperature  law  are  equiva¬ 
lent,  and  that  they  are  simply  related  through  the  concept  of  a  ref¬ 
erence  temperature. 
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VELOCITY  AND  ENTHALPY  PROFILES  WITH  A  POWER-LAW  VISCOSITY 

The  effects  of  changing  the  exponent  in  the  power- law  expression 
for  the  viscosity  are  shown  in  Figs.  1  and  2.  A  value  uo  *  0.7  cor¬ 
responds  closely  to  conventional  wind-tunnel  conditions,  while  uu  *  0.5 
represents  conditions  encountered  in  hypersonic  flight.  The  value 
a)  *  1  represents  the  simplification  (pp.)  «  constant  which  has  been 
widely  used  in  studying  the  laminar  boundary  layer.  The  values 
0  *  0.5  and  a  =  0  are  appropriate  to  axisymmetric  stagnation-point 
flow,  while  0  =  0  and  a  ■  1  represent  the  hypersonic  laminar  boundary 
layer  over  a  flat  plate. 

Several  important  effects  may  be  discerned  in  Figs.  1  and  2.  First, 

★ 

decreasing  values  of  u)  decrease  the  0.90  thicknesses  of  the  velocity 
and  total-enthalpy  boundary  layers.  Second,  the  variable  viscosity 
introduces  a  distinct  inflection  in  the  velocity  and  enthalpy  profiles, 
as  evidenced  by  a  maximum  in  the  derivatives  f '  (T[)  and  0  *  01)  .  Since 
the  existence  of  an  inflection  point  in  the  velocity  profile  is  re¬ 
lated  to  the  stability  of  the  laminar  boundary  layer,  this  effect  may 
explain  the  anomalous  result  of  stability  theories  based  on  uu  *  1  that 
the  critical  transition  Reynolds  number  for  a  flat  plate  becomes  infi¬ 
nite  with  a  highly  cooled  wall.  The  fact  that  the  wall  derivatives 
and  0^  for  0*0  are  nearly  independent  of  the  value  of  u#  appears  to 
be  fortuitous.  The  variations  of  f 1  ,  f"  ,  0,  and  9*  with  uu  are  consid¬ 
erably  reduced  at  higher  values  of  the  wall  temperature  t^. 


EFFECTS  OF  VISCOSITY  LAW  ON  SKIN  FRICTION  AND  HEAT  TRANSFER 

One  of  the  purposes  of  this  study  is  to  examine  the  use  of  the 

Sutherland  viscosity  law  in  determining  the  properties  of  the  hypersonic 

laminar  boundary  layer.  In  Fig.  3,  the  effects  of  pressure  gradient, 

wall  temperature,  and  the  Sutherland  constant  8  on  the  normalized  heat 

transfer  9'  are  shown.  The  normalized  heat  transfer  9'  is  defined  by 
w  w 


<1  - 

(taw  ‘  SP 


(10) 


*The  points  where  the  velocity  and  total  enthalpy  reach  0.90  of 
their  values  outside  the  boundary  layer. 


(a)  Velocity  and  shear  distributions 


(b)  Entholpy  ond  enthalpy-gradient  distributions 


Fig.  1 — Velocity  and  enthalpy  distributions  across  the  laminar 
boundary  layer:  /9  ■  0,<7- 1 


V 

(b)  Enthalpy  and  enthalpy-gradient  distributions 


Fig.  2— -Velocity  and  enthalpy  distributions  across  the  laminar 
boundary  layer:  ■  0. 5,  <r  •  0 
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where 

t  =  adiabatic  wall  temperature 
t^  »  wall  temperature 

0^  =  total  enthalpy  gradient  at  the  wall. 

The  Blasius  result  (0=0,  Pr  *  1,  pp.  ■  constant)  is  given  for  compari¬ 
son.  The  Sutherland  constant  s  is  defined  by 


o 


where  S  is  a  reference  temperature  for  the  particular  gas  under  consid¬ 
eration.  The  constant  s  **  0.2  for  ordinary  wind-tunnel  conditions,  and 
s  *  0.02  for  hypersonic  flight  conditions. 

The  pressure-gradient  effect  is  as  expected.  If  the  wall  tempera¬ 
ture  and  Sutherland  constant  are  the  same,  then  an  increase  of  the  pres¬ 
sure  gradient  increases  the  heat  transfer  at  the  wall.  The  dependence 
of  the  heat  transfer  on  the  Sutherland  constant  increases  as  the  wall 
temperature  decreases.  At  a  value  of  t^  *  0.6,  the  normalized  heat 
transfer  0^  is  essentially  independent  of  s.  If  the  wall  temperature 
is  reduced  to  t^  -  0.15,  then  0^  depends  strongly  on  s;  the  normalized 
heat  transfer  3^  increases  with  increasing  s. 

Figure  4  shows  the  variation  of  the  normalized  heat  transfer  0^ 
as  a  function  of  the  power- law  exponent  u>  for  the  hypersonic  limit 
o  m  l.  Here  again,  an  increase  in  the  pressure  gradient  or  the  wall 
temperature  results  in  an  increase  in  the  normalized  heat  transfer. 

The  values  of  0^  seem  to  be  more  sensitive  to  the  pressure  gradient  in 
the  hypersonic  limit  a  *  1  and,  in  addition,  are  somewhat  higher  than 
they  are  at  a  *  0.5.  In  general,  the  maximum  variation  of  0^  at  a 
given  wall  temperature  and  pressure  gradient  with  u)  or  s  is  less  than 
30  per  cent. 

Figures  5  and  6  are  included  to  illustrate  the  magnitude  of  the 

effects  of  the  viscosity- temperature  relation  on  skin  friction  and 

heat  transfer.  In  Fig.  5,  the  normalized  heat  transfer  O'  is  shown 

w 

as  a  function  of  t^  and  uu.  It  is  interesting  to  note  that  the  single 


0.7 


Pr  =  0.7 


0-3  h 

0 1 _ I _ I _ I _ I - 1 - 

0  0.2  0.4  0.6  0.8  1.0  1.2 

Wall  temperature,  t„ 

Fig.  5 — Effects  of  wall  temperature  on  heat  transfer 


0  0.1  0.2  0.3  0.4  0.5  0.6 

Pressure-gradient  parameter,  )9 

Fig.  6 — Effects  of  pressure  gradient  on  skin  friction 
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Blasius  value  0^  *  0.4696  for  Pr  ■  u)  =  1  and  P  *  0  is  a  fortuitously 
good  approximation  to  the  exact  value  of  0^  over  a  large  range  of 
conditions.  For  any  particular  choice  of  conditions,  however,  such  a 
simple  approximation  may  be  in  error  by  as  much  as  30  per  cent.  The 
changes  in  0^  introduced  by  variation  of  the  exponent  u)  are  of  the  same 
order  as  those  introduced  by  the  pres sure- gradient  parameter  P  and  the 
wall  temperature  t^.  It  is  important,  therefore,  that  a  proper  vis¬ 
cosity  formulation  be  used  when  exact  solutions  are  desired. 

The  wall  velocity  gradient  f”  is  strongly  affected  by  pressure 
gradient,  as  shown  in  Fig.  6.  This  is  because  the  parameter  P  appears 
explicitly  in  the  momentum  equation  (1) ;  the  effects  of  pressure  grad¬ 
ient  appear  in  the  energy  equation  (2)  only  indirectly  through  the  ve¬ 
locity  ratio  f’  *  (u/ug) .  It  may  be  seen  in  Fig.  6  that  the  effects 
of  the  viscosity  exponent  on  f^  are  small  compared  to  the  effects  in¬ 
troduced  by  the  pres sure- gradient  parameter.  The  effects  of  large  val¬ 
ues  of  P  on  the  wall  shear  f"  have  been  studied  extensively  by  Cohen 

(4)  w 

and  Reshotko'  y  and  others. 

RELATION  BETWEEN  SUTHERLAND  AND  POWER-LAW  VISCOSITY  FUNCTIONS 

Another  purpose  of  the  Sutherland- law  investigation  was  to  deter¬ 
mine  the  relations  between  solutions  obtained  using  a  Sutherland- law 

U) 

viscosity  and  those  obtained  by  the  use  of  a  power- law  ^  ~  T  employing 
an  empirically  calculated  u^.*  Solutions  for  specific  boundary  condi¬ 
tions  were  obtained  using  the  Sutherland- law  viscosity  relationship. 

A  reference  temperature,  t^,  was  then  used  which  was  originally  sug¬ 
gested  by  Rubesin  and  Johnson^  and  modified  by  Eckert^  to  provide 
an  approximation  method  of  evaluating  the  wall  conditions  for  a  compres¬ 
sible  boundary  layer.  The  Eckert  reference  temperature  is  given  by 

tr  -  0.5  <te  +  t^  +  0.22  (1  -  y  (12) 

*It  should  be  emphasized  that  the  Sutherland  law  itself  is  de¬ 
rived  from  an  empirical  approximation  to  the  intermolecular  potential. 
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and  was  used  to  calculate  an  empirical  u>^  from  the  following  formula: 


An  exact  solution  to  the  similar  equations  was  then  obtained  using  the 

^r 

temperature-viscosity  relation  \x  ~  T  ,  and  these  results  were  compared 
with  solutions  obtained  using  the  Sutherland  law.  Figure  7  shows  the 
differences  in  the  wall  derivative  f^  obtained  for  some  typical  cases. 
Specifically,  the  graph  shows  Af^  *  f^(s)  "  (oo  )  plotted  as  a  func¬ 
tion  of  s  for  different  wall  temperatures  and  pressure  gradients.  The 
maximum  error  is  approximately  2.5  per  cent  in  the  wall  velocity  grad¬ 
ient;  in  most  cates,  the  error  is  well  below  1  per  cent.  Figure  7  also 
shows  the  effects  of  the  important  parameters  t  and  P  on  the  differ¬ 
ences  Af1^  induced  by  the  use  of  the  calculation.  Clearly,  the  dif¬ 
ferences  increase  as  s  increases,  and  both  the  difference  function  and 
s  must  go  to  zero  together  is  Identically  1/2  for  s  *  0)  .  For  the 
cold  wall,  the  difference  is  greater  because  the  power- law  approxima¬ 
tion  cannot  fit  the  Sutherland -law  viscosity  throughout  the  larger  tem¬ 
perature  range.  These  differences  are  reduced  to  negligible  values 
for  t^  ^  0.6.  The  effect  of  pressure  gradient  at  low  wall  temperature 
is  also  evident:  An  increase  in  P  gives  rise  to  an  increase  in  the 
difference  function  Af^.  The  higher  pressure  gradient  causes  the  power- 
law  approximation  to  be  a  less  accurate  fit  of  the  Sutherland  law.  In 
general,  the  differences  are  quite  small  even  for  the  cold-wall  cases, 
and  the  power-law  approximation  using  the  Eckert  reference  temperature 
should  prove  to  be  a  valuable  empirical  method. 

The  values  of  u)  as  a  function  of  the  Sutherland  constant  8  are 
r 

shown  in  Fig.  8  for  the  three  wall  temperatures  t  *  0.15,  t^  -  0.4, 

and  t  *0.6. 

w 

It  can  be  seen  that  the  empirically  determined  u»  is  a  smooth 
function  of  the  Sutherland  constant  s,  so  it  will  be  possible  to  in¬ 
terpolate  for  conditions  not  presented  in  this  Memorandum.  When 


Fig.  7 — Wall  velocity-gradient  differences  between 
the  two  viscosity  laws 
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IV,  EFFECTS  OF  THE  PRANDTL  NUMBER 


The  effects  of  the  Prandtl  number  on  the  enthalpy  and  velocity 
distributions  in  the  boundary  layer  are  shown  in  Fig.  9  for  the  cold- 
wall  case  (t^  =  0.15)  in  the  hypersonic  limit,  a  *  1.  The  velocity 
profiles  are  not  noticeably  affected,  but  there  is  a  drastic  change  in 
the  enthalpy  gradient  as  well  as  the  enthalpy  distribution  within  the 
boundary  layer.  The  wall  enthalpy  gradient  is  considerably  reduced 
when  the  Prandtl  number  is  decreased  from  1.0  to  0.5.  The  enthalpy 
profile  is  less  full  for  the  Pr  *  0.5  case,  and  this  has  the  effect  of 
concentrating  the  high  enthalpy  gradients  near  the  outer  edge  of  the 
boundary  layer. 

In  order  to  demonstrate  the  effects  of  the  Prandtl  number  more 

decisively,  the  wall  enthalpy  gradient  0^  and  the  normalized  heat- 

transfer  parameter  are  plotted  as  functions  of  Pr  in  Fig.  10.  It 

is  shown  in  (a)  of  Fig.  10  that  0^  increases  with  increasing  Pr;  here, 

a  ®  0,  so  that  the  dissipation  term  appearing  in  the  energy  equation 

(2)  is  identically  zero.  It  may  be  seen  that  the  common  assumption  of 

Pr  «  1  for  air  introduces  an  error  in  0*  which  is  of  the  same  order  of 

w 

magnitude  (about  5-15  per  cent)  as  the  change  in  0^  due  to  a  varia¬ 
tion  of  the  pressure-gradient  parameter  0  from  0  to  0.5.  It  is  there¬ 
fore  numerically  inconsistent  (at  least  for  o  *  0)  to  include  the  ef¬ 
fects  of  pressure  gradient  while  neglecting  the  effects  of  Prandtl  num¬ 
ber  . 

In  the  hypersonic  limit,  o  ■  1,  the  situation  changes  dramatically; 
(b)  of  Fig.  10  indicates  the  variation  of  0^  with  Prandtl  number  for 
P  -  0  with  and  without  the  viscous-dissipation  term.  The  effects  of 
Pr  and  wall  temperature,  t^,  on  01  are  greatly  Increased  by  including 
the  viscous-dissipation  term.  This  is  primarily  because  of  the  large 
change  in  the  adiabatic  wall  temperature  t ^  with  Prandtl  number.  A 
useful  empirical  formula  which  relates  the  adiabatic  wall  temperature 
t^w  to  the  parameters  Pr,  u),  a,  g,  and  f  is 

t  -  1  +  1.428  o  (,/Pr  -  1)  (1  +  0.290)  (1  -  0.3u>)  (1  -  0.95^. 

(14) 


(o)  Velocity  and  velocity-gradient  distributions 


V 

(b)  Enthalpy  and  enthalpy-gradient  distributions 


Fig.  9 — Velocity  and  enthalpy  distributions  for  different 
Prandtl  numbers 
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Prondtl  number,  Pr 


Prandtl  number,  Pr 


Prandtl  number,  Pr 

Fig.  10 — Effects  of  Prandtl  number  on  heat  transfer 


o* 


-19- 


Equation  (14)  becomes  somewhat  inaccurate  for  values  of  f  :  <  -0.6  and 
3  >  0.5,  but  for  moderate  values  of  these  parameters  the  predicted  val¬ 
ue  of  t  is  correct  to  better  than  1  per  cent.  From  Eq.  (14)  and  (b) 
of  Fig.  10,  the  important  effects  of  Prandtl  number  on  surface  heat 
transfer  may  be  seen.  The  familiar  Blasius  value  0^  ■  0.4696  for  a 
flat  plate  with  Pr  =  uu  =  1  may  overestimate  the  surface  heat  transfer 
in  hypersonic  flow  by  a  factor  of  2  or  more,  even  for  a  moderately 
cooled  wall . 

The  normalized  enthalpy  gradient  0'  shows  a  much  smaller  change 

w 

with  Prandtl  number  than  the  gradient  0'.  It  is  indicated  in  (c)  of 

w 

Fig.  10  that  the  normalized  derivative  0'  is  close  to  the  classical 
°  —  w 

Blasius  value  even  in  hypersonic  flow.  From  Figs.  3,  4,  5,  and  10, 
it  may  be  concluded  that  classical  value  0'  =  0.4696  is  a  good  approx- 
imation  over  large  ranges  of  Prandtl  number,  viscosity  parameter,  pres¬ 
sure-gradient  parameter,  and  wall  temperature.  Combining  this  result 
with  Eq.  (14),  the  enthalpy  derivative  0^  may  be  determined. 
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v.  effects  of  large  injection 


Early  theoretical  calculations  for  mass  transfer  into  a  laminar 

(7  8) 

boundary  layer  on  a  flat  plate'  *  7  indicated  the  existence  of  "blow- 
off"  points,  i.e.,  separation  of  the  boundary  layer,  at  high  mass-in- 
jectiijn^rates .  These  points  were  determined  by  noting  the  values  of 

fw  ■  cTu2  y  (Rex)  for  which  £w  ~  (^y),  *  °*  Early  effort8  t0  cor- 

relateCt^e  exact  solutions  for  the  laminar  boundary  layer  on  a  flat 

(9) 

plate  with  mass  transfer  (see,  e.g.,  Gross,  et  al . , v  7)  suggested  the 
use  of  a  linear  relationship  between  the  mass-in jection  parameters  f^ 
and  the  dimensionless  skin  friction: 


(15) 


where 

c^  ■  local  skin  friction  with  mass  transfer 

c^  «  local  skin  friction  without  mass  transfer, 
o 

Extrapolation  of  this  result  would  give  an  f w  «  -  0.5  for  the  blowoff 

point  where  cf  -  0.  Hartnett  and  Eckert,  Scott,  ^  and  Hayday^^ 

showed  that  a  similar  blowoff  point  was  not  evident  in  the  case  of 

(9) 

stagnation-point  flows.  In  Gross,  et  al .  , v  7  a  comparison  is  made  be¬ 
tween  flat-plate,  wedge,  and  stagnation-point  flows  to  show  the  effects 
of  pressure  gradient  on  the  dimensionless  heat  transfer  0^.  The  im¬ 
plication  was  that  a  blowoff  point  for  stagnation  flows--if  it  existed-- 
would  occur  at  a  much  greater  value  of  f^  than  for  the  flat  plate. 

In  this  study,  calculations  were  carried  out  to  determine  the  be¬ 
havior  of  the  velocity  gradient  and  the  enthalpy  gradient  at  the  wall 
for  a  stagnation-point  boundary  layer  for  large  values  of  the  blowing 
parameter  |  |.  The  results  shown  in  Fig.  11  indicate  that  the  solu- 

tions  are  asymptotic  and  that  a  blowoff  point  is  never  achieved,  even 
though  the  values  of  the  enthalpy  gradient  become  quite  small.  It 
seems  clear  from  the  calculations  that  will  approach  zero  only  at 
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"f* 

Fig.  11 — Wall  velocity  gradient  and  heat  transfer  with 
large  injection:  two-dimensional  stagnation  point 
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very  large  values  of  -f^.  Figures  12  and  13  show  the  development  of 
the  profiles  for  the  velocity  and  enthalpy  and  their  gradients  for  in¬ 
creasing  values  of  the  mass- transfer  parameter  |  f^  |.  The  solid-wall 
profiles  have  been  included  for  comparison.  It  is  interesting  to  note 
the  drastic  effect  of  the  mass- transfer  parameter  in  shifting  the  peak 
enthalpy  toward  the  outer  edge  of  the  boundary  layer.  The  increase  in 
the  thickness  of  the  boundary  layer  with  mass  injection  is  graphically 
illustrated  in  these  figures. 
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Fig.  12 — Velocity  and  velocity-gradient  distributions 
with  injection 
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Fig.  13 — Enthalpy  and  enthalpy-gradient  distributions 
with  injection 
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VI.  BEHAVIOR  OF  THE  BOUNDARY  LAYER  IN  THE  HYPERSONIC  LIMIT 

2  2 

The  influence  of  the  hypersonic  parameter  o  *  (u^  /2He)  (ue/uoo) 
on  the  wall  shear  f^,  the  normalized  enthalpy  gradient  0^,  and  the  dis¬ 
placement— thickness  integral  1^  is  considered  in  Figs.  14  and  15.  The 
ordinate  is  given  in  terms  of  the  ratios 


f" 

w 


1 


1 


where  the  subscript  a  =  1  indicates  a  value  for  the  hypersonic  limit 

a-*l[T/T)->0].  The  integral  I.  is  related  to  the  displacement 

e  $  (1)  1 
thickness  6  by  the  equations  ' 


6 


* 


lili 

PwUe 


1/2 


h  *l[(1  '  f’2>  -  d  -  V(1  '  0)]  dTlj 

O 


(16) 

(17) 


I2  -  J  f’(l  -  f')  dH.  (18) 

o 

For  the  hypersonic  boundary  layer,  (TQ/Te)  ^  2  ")  >:>  ^ 9  aru*  ^ 

is  very  nearly  proportional  to  the  product  Me  1^ . 

The  ratios  of  f^,  0^,  and  1^  increase  with  increasing  values  of 

a.  The  maximum  change  in  the  ratios  is  about  20  per  cent  between  o  -  0 

and  a  »  1.  One  important  conclusion  may  be  inmediately  drawn:  The 

2 

similarity  requirement  that  the  ratio  (u^/u^)  be  constant  along  the 
body  surface  is  not  a  severe  restriction,  since  large  changes  in  this 
quantity  (i.e.,  in  a)  do  not  have  a  proportionately  large  effect  on 
the  boundary- layer  solutions. 
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Fig.  14 — Influence  of  hypersonic  parameter  a  on  skin  friction, 
heat  transfer,  and  displacement  thickness:  tw  •  0. 15 


Fig.  15—  Influence  of  hypersonic  parameter  a  on  skin  friction, 
heat  transfer,  and  displacement  thickness:  tw  -  0.60 
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The  change  in  9^,  and  1^  with  o  is  nearly  independent  of  P. 

A  comparison  of  Figs.  14  and  13  shows  that  the  ratios  suffer  a  larger 
change  with  a  lower  wall  temperature.  From  Eq.  (14),  and  Figs.  14 
and  15,  it  may  be  deduced  that  the  enthalpy  derivative  9^  decreases 
with  increasing  o,  and  the  variation  is  greater  with  increasing  wall 
temperature. 
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VII.  CONCLUDING  REMARKS 


The  results  of  this  study  may  be  succinctly  stated  as  follows: 

1.  The  effect  of  changing  u>  in  the  power-law  viscosity  p.  a  T^ 
on  the  enthalpy  and  velocity  profiles  is  most  important  for  low  wall 
temperatures.  Differences  of  20  -  30  per  cent  in  predicted  heat  trans¬ 
fer  and  skin  friction  for  0.5  ^  u)  £  1  are  observed  at  t  »  0.15. 

w 

2.  The  effect  of  varying  the  Prandtl  number  on  the  enthalpy  dis¬ 
tribution  in  the  boundary  layer  is  important  for  low  wall  temperature 
(t^  ^  0.15)  and  large  values  of  the  hypersonic  parameter  a. 

3.  The  laminar  boundary  layer  at  a  stagnation  point  for  a  cold 
wall  does  not  "blow  off"  or  separate  even  at  high  blowing  rates.  Cal¬ 
culations  indicate  that  0^  -*  0  much  more  rapidly  than  f^,  and  that  for 

f  <  -  1.2,  0'  -  0,  while  f"  is  finite, 
w  w  w 

4.  The  numerical  solutions  show  that  the  viscous-dissipation 
* 

term  has,  at  most,  a  10  -  20  per  cent  effect  on  the  wall  shear  f^, 
the  displacement  thickness  1^,  and  the  normalized  enthalpy  gradient 
9^;  these  three  quantities  reach  a  maximum  value  for  a  -•  1,  i.e.,  in 
the  hypersonic  limit. 

5.  If  the  enthalpy  gradient  0^  is  normalized  to  0^  where 


0 1 
w 


» -  v 

<‘.v  -  V 


then  the  function  9^  is  only  weakly  dependent  upon  the  Sutherland  con¬ 
stant  s.  In  general,  0^  increases  with  both  wall  temperature  t^  and 
pressure-gradient  parameter  0. 

6.  The  boundary- layer  characteristics  calculated  using  a  Suther- 
land-law  viscosity  can  be  approximated  almost  exactly  using  a  pover- 
law  viscosity  p  a  T  if  the  is  determined  empirically  from  the 
equation 


Third  term  on  left  side  of  Eq.  (2). 
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where  is  the  Eckert  reference  temperature,  and  b  is  a  temperature 
constant  for  the  gas. 
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